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Abstract 

Deep-inelastic positron-proton interactions at low values of Bjorken-x down to x rs 
4 • 10~ 5 which give rise to high transverse momentum 7r°-mesons are studied with the 
HI experiment at HERA. The inclusive cross section for ir° -mesons produced at small 
angles with respect to the proton remnant (the forward region) is presented as a function 
of the transverse momentum and energy of the ir° and of the four-momentum transfer Q 2 
and Bjorken-x. Measurements are also presented of the transverse energy flow in events 
containing a forward 7r°-meson. Hadronic final state calculations based on QCD models 
implementing different partem evolution schemes are confronted with the data. 
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1 Introduction 



Measurements of the hadronic final state in deep-inelastic positron-proton scattering (DIS) at 
HERA have allowed precision tests of the theory of the strong force, Quantum Chromodynam- 
ics (QCD). In particular, properties of the hadronic final state in the region near to the proton 
remnant system (hereafter referred to as the 'forward region') have been shown to be sensitive 
to the QCD radiation pattern formed from the cascade initiated by a parton from the proton be- 
fore it undergoes a hard scatter [1-4]. This paper presents studies made by the HI experiment 
of DIS interactions at values of Bjorken-rr down to x ~ 4 • 10~ 5 containing at least one high 
transverse momentum, forward going 7r°-meson. 

A generic diagram for parton evolution in a DIS process at low x in which a gluon from 
the proton undergoes a QCD cascade is shown in Fig. 1 (a). The gluon eventually inter- 
acts with the virtual photon via a hard photon-gluon fusion process which can be calculated 
within perturbative QCD using an exact matrix element. Several perturbative QCD-based 
prescriptions are available to describe the dynamics of the parton evolution process. The 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [5] resum leading 
\og(Q 2 ) terms and ignore log(l/x) terms. In an axial gauge this corresponds to the resum- 
mation of diagrams in which the parton cascades follow a strong ordering in transverse mo- 
menta k\ n 3> k\ n _ x 3> • • ■ 3> k\ x . At sufficiently small values of x the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation [6] should be applicable since the log(l/x) terms should 
dominate the evolution. In this scheme the cascade is ordered strongly in fractional momenta 
x n <C x n -\ C ••• < xi, while the transverse momenta perform a 'random walk' with kn 
being close to kn-i, though it can be both larger or smaller. The Ciafaloni-Catani-Fiorani- 
Marchesini (CCFM) equation [7] interpolates between the DGLAP and BFKL approximations 
with parton emissions ordered in angle. 

A parton chain without a requirement of kr ordering throughout the complete cascade is 
provided in a picture of low x DIS in which the virtual photon is ascribed a partonic structure. 
This process is illustrated in Fig. 1 (b) in which /c T -ordered DGLAP cascades are initiated both 
from the proton and photon, leading to the hard interaction at the centre of the QCD 'ladder' . 

Hadronic final state observables are sensitive to the dynamics of QCD processes and are thus 
expected to be able to discriminate between different evolution approximations. As illustrated 
in Fig. 1, the selection of leading particles or jets in the forward region can tag individual 
parton emissions at high transverse momentum. An advantage of using single particles rather 
than jets is that ambiguities due to the choice of jet algorithm are removed. On the other hand, 
uncertainties due to hadronisation are typically larger for single particle studies. As was pointed 
out in [8], the selection of particles or jets with values of transverse momentum squared of 
similar magnitude to Q 2 suppresses the contribution of /c r -ordered cascades with respect to k T - 
unordered processes. In addition, the phase space for BFKL effects is enhanced if the fraction 
of the proton's energy of the particle or jet is required to be greater than Bjorken-x. 

A recent HI study of forward going 7r°-mesons [9] found that a model which implemented 
a DGLAP parton cascade from the proton significantly underestimated the cross section at low 
values of Bjorken-rr. Leading order BFKL calculations with kinematic constraints which mimic 
higher orders and a model implementing virtual photon structure gave a better description of 
the data. 
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The pseudorapidity dependence of the mean transverse energy produced in an event (the so- 
called transverse energy flow) provides a complementary means of studying the hadronic final 
state. Compared with studies of jets and high transverse momentum particles, measurements 
of transverse energy flow typically cover a wider range of pseudorapidity and are sensitive to 
parton emissions of lower transverse momentum. Previous measurements of transverse energy 
flow at HERA were indeed found to be sensitive to the modelling of both the perturbative QCD 
evolution and the soft hadronisation process [1,10, 11]. Measurements of transverse energy flow 
in events containing particles with high transverse momentum also reveal the range over which 
transverse momentum is compensated following the emission of QCD radiation [12]. 

This paper presents a study of low x DIS interactions in which high transverse momentum 
7r°-mesons are produced in the forward region. The results are based on a data sample which is 
more than three times larger than that used for earlier studies [9]. Consequently, the inclusive 
7r° cross section is measured with greater precision and more differentially as a function of x, 
Q 2 and the transverse momentum and energy of the 7r°-meson. Furthermore, for the first time, 
measurements are presented of transverse energy flow for ep interactions containing forward 
going 7r°-mesons. This allows a more complete investigation of hadronic final states containing 
a hard forward ir° than was previously possible. 



2 QCD-based Models 

Calculations of the expected production rate of 7r°-mesons and the associated transverse energy 
flow in DIS are available in the form of Monte Carlo event generators. These use first-order 
QCD matrix elements and adopt various approaches to modelling the parton cascade. These 
models are used to provide comparisons of theory predictions with the measurements presented 
here. In addition they are used, together with a simulation of the HI detector, to correct the 
measurements for the finite acceptance and resolution of the detector. Unless otherwise stated 
the proton and virtual photon parton densities used in these models are CTEQ6M [13] and 
SAS-1D [14], respectively. 

Lepto 6.51 [15] matches first-order QCD matrix elements to DGLAP-based leading-log 
parton showers. The factorisation and renormalisation scales are set to Q 2 . Lepto also al- 
lows for non-perturbative rearrangement of the event colour topology via so-called soft color 
interactions in the final state [16]. 

Rapgap 2.08/20 [17] also matches first-order QCD matrix elements for direct photon pro- 
cesses to DGLAP-based leading-log parton showers. In addition to the direct photon processes, 
Rapgap simulates resolved photon interactions in which the virtual photon is assumed to have 
partonic structure. For the predictions presented here, the renormalisation and factorisation 
scales are set to + Q 2 , where px is the transverse momentum of the partons emerging from 
the hard scattering process. The hadronic final state predictions of Rapgap, when only direct 
photon interactions are considered, are very similar to those of Lepto. 

Ariadne 4. 10 [18] is an implementation of the Colour Dipole Model (CDM) [19] of a chain 
of independently radiating dipoles formed by emitted gluons. Since all radiation is assumed to 
come from the dipole formed by the struck quark and the remnant, photon-gluon fusion events 
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have to be added and are taken from the QCD matrix elements. The parameters in Ariadne 
have been optimised in order to describe a range of hadronic final state measurements [20]. 

In order to study the effects of initial and final state QED radiation, the above models are 
interfaced with Heracles [21] within the Django [22] model. 

Cascade 1.0 [23] uses off-shell QCD matrix elements, supplemented with parton emis- 
sions based on the CCFM equation within a backward evolution approach. An unintegrated 
gluon density, obtained using CCFM evolution and fitted to describe the inclusive DIS cross 
section [24], is used as an input to this model. In the present analysis an updated version [25] 
of Cascade with an improved treatment of the soft region and a new parameterisation of the 
unintegrated gluon density is used. These modifications provide an improved description of 
forward jet production [26]. 

To perform the hadronisation step, all of the above models use the LUND string fragmenta- 
tion [27] scheme, as implemented in Jetset [28] in case of Lepto, Rapgap and Ariadne 
and in Pythia [28] for Cascade. 

Predictions of the n° cross sections are also available from an analytical calculation at the 
parton level [29] based on a modified BFKL evolution equation at lowest order. These calcula- 
tions are then convoluted with a n° fragmentation function [30]. The proton parton densities are 
taken from [31]. The modified evolution equation imposes a "consistency constraint" [32, 33] 
which, it is argued, mimics much of the contribution from non-leading log(l/a;) terms to the 
BFKL equation. However, these predictions are very sensitive to the choice of scale for the 
strong coupling constant a s and to the infra-red cut-off. In the present analysis, the scale for 
a s is taken to be the squared transverse momentum of the emitted partons, k\, and the infrared 
cut-off in the modified BFKL equation is set at 0.5 GeV 2 . 

Recently, calculations of the cross section for the production of high transverse momen- 
tum hadrons in DIS interactions have been made [34], which describe earlier measurements of 
forward ir° production [9]. These comprise next-to-leading order (NLO) matrix elements, con- 
voluted with NLO fragmentation functions [30]. The proton parton densities MRST99 (higher 
gluon) [35] are used. For comparison with the HI data, the renormalisation, factorisation and 
fragmentation scales are each set to (Q 2 + p^>)/2, where is the transverse momentum, in 
the photon-proton centre of mass system 1 , of the parton which fragments into the forward ir°. 
In these calculations, a large part of the cross section is generated by higher order contributions 
which correspond to lowest order BFKL and resolved photon processes [34]. 



3 Experimental Apparatus 

A detailed description of the HI detector can be found elsewhere [36]. The following section 
briefly describes the components of the detector which are most relevant for this analysis. 

A liquid argon (LAr) calorimeter is used to measure the hadronic energy flow and the can- 
didate 7r°-meson properties. The LAr calorimeter provides measurements over the laboratory 

'All quantities presented in the hadronic centre-of-mass system are denoted by the superscript *. 
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polar angle range 4° < 6 < 154°, where 9 is defined with respect to the direction of the proton 
beam, and offers full azimuthal coverage. It consists of an electromagnetic section with lead 
absorbers and a hadronic section with steel absorbers. Both sections are highly segmented in 
the transverse and longitudinal directions. The total depth of both sections varies between 4.5 
and 8 interaction lengths in the region 4° < 9 < 128°, and between 20 and 30 radiation lengths 
in the region 4° < 9 < 154° increasing towards the forward direction. The fine granularity of 
the electromagnetic section in the forward direction is characterized by four-fold longitudinal 
segmentation and a typical lateral cell size of 3.5 x 3.5 cm 2 . Test beam measurements of the 
LAr calorimeter modules showed an energy resolution of a E /E pa 0.50/ \/E[GeV] © 0.02 for 
charged pions and of cr E /E pa 0.12/y/E[GeV\ © 0.01 for electrons [36]. The hadronic energy 
measurement is made by applying a weighting technique to the electromagnetic and hadronic 
components of the energy deposition, in order to account for the non-compensating nature of 
the calorimeter. The absolute scales of hadronic and forward-reconstructed electromagnetic 
energies are known to 4% [37] and 3% [38], respectively. 

The SPACAL is a lead/scintillating fibre calorimeter covering the region 153° < 9 < 177.5° 
with an electromagnetic and a hadronic section. It is used to measure the scattered positron en- 
ergy and the backward hadronic energy flow. The energy resolutions for electrons and hadrons 
are a E /E pa 0.07/ y/E[GeV] © 0.01 [39] and a E /E pa 0.3/ ^/E [GeV] [40], respectively. The 
energy scale uncertainties are 1% for electrons and 7% for hadrons [41]. 

The calorimeters are surrounded by a superconducting solenoid which provides a uniform 
magnetic field of 1.15 T in a direction parallel to the proton beam in the tracking region. 
Charged particle tracks are measured in the central tracker (CT) and forward tracker (FT) sys- 
tems which cover the polar angle ranges of 25° < 9 < 155° and 5° < 9 < 25°, respectively. 
Information from the CT is used in this work to trigger events, to locate the event vertex and 
also contributes to the measurement of transverse energy. 

A backward drift chamber (BDC) in front of the SPACAL with an angular acceptance of 
151° < 9 < 177.5° serves to identify electron candidates and to precisely measure their di- 
rection. Using information from the BDC, the SPACAL and the reconstructed event vertex 
position, the polar angle of the scattered electron is known to about 0.5 mrad [41]. 

The luminosity is measured using the Bethe-Heitler process ep — > ejry with two TICl/TlBr 
crystal calorimeters installed in the HERA tunnel. 

4 Data Analysis 

The data used for this analysis were collected in 1996 and 1997 when positrons and protons 
with energies of 27.6 GeV and 820 GeV, respectively, were collided. The data-set used in this 
work corresponds to an integrated luminosity of 21.2 pb _1 . 

4.1 DIS Event Selection 

DIS events are selected by triggers based on electromagnetic energy deposits in the SPACAL 
calorimeter and the presence of charged particle tracks in the CT. For the 7r°-enriched event 
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sample, the trigger efficiency lies between 60% and 80%, determined using independently trig- 
gered data. The inefficiency is mainly due to the suppression by the trigger of events with less 
than 3 charged particles measured in the CT. The data are corrected for the trigger inefficiency 
by applying a weight to every selected DIS event. 

The event kinematics are calculated from the polar angle and the energy of the scattered 
positron. In order to maintain optimal efficiency and acceptance, scattered positron candidates 
in the SPACAL are required to have an energy E e > > 10 GeV and to lie in the region of 
polar angle of 156° < 9 e > < 177°. The data-set is further restricted to the kinematic range in 
inelasticity y and virtuality Q 2 of 0.1 < y < 0.6 and 2 < Q 2 < 70 GeV 2 , respectively. The 
resulting values of Bjorken-rr extend over two orders of magnitude in the range 4 ■ 10~ 5 < x < 
6 • 10~ 3 . 

To further reject background from photoproduction interactions in which a particle from 
the hadronic final state is misidentified as a scattered positron in the SPACAL, the condition 
35 < T,j{Ej — p z j) < 70 GeV is applied. Here Ej and p z j are the energy and longitudinal 
momentum, respectively, of a particle, and the sum extends over all particles in the event except 
those detected in the luminosity system. 

4.2 Forward 7r° Selection 

The selection of forward 7r°-mesons closely follows that of earlier work [3, 9,42]. The tt° can- 
didates are identified via the dominant decay channel ir° — > 27 using calorimetric information 
only. Criteria are placed on the kinematic properties of the candidate to ensure high acceptance 
and efficient background rejection. The energy E n of the n° scaled with the energy E p of the 
proton beam, x w = E n / E p , is required to be greater than 0.01 . As a consequence of this cut, the 
decay photons are not resolved individually but are merged into a single electromagnetic cluster 
in the detector. The candidates are also required to lie in the polar angle region 5° < 9 W < 25°. 
This polar angle region, referred to as 'forward' in the laboratory frame, corresponds to the 
central region in the hadronic centre of mass frame, —1.25 < 77* < 2.0. Furthermore, the n° 
transverse momentum in the photon-proton centre of mass system, p^^ must be greater than 
2.5 GeV. The Lorentz boost to this frame is calculated using the kinematics of the scattered 
electron. In this frame the proton direction is chosen to define the negative z* axis. 

Electromagnetic and hadronic showers are discriminated by an analysis of the longitudinal 
and transverse shapes of the energy depositions. A ir° candidate is required to have more than 
90% of its energy deposited in the electromagnetic part of the LAr calorimeter. A "hot" core 
of the most energetic group of contiguous cells within a cluster, which must include the hottest 
cell, is required to account for over 50% of the cluster energy. The lateral spread of the shower, 
defined as in [43], is required to be less than 4 cm. Electromagnetic clusters with a small 
longitudinal extent are selected by requiring that the difference in the amounts of energy found 
in the second and fourth longitudinal layers of the electromagnetic calorimeter must be more 
than 40% of the total cluster energy. 

Following this selection approximately 5500 (2000) ir° candidates remain after the trans- 
verse momentum cut p^ w > 2.5 (3.5) GeV. Using simulated events generated with the Lepto 
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and Ariadne models the efficiency of the selection is estimated to be approximately 45%, and 
the contribution of non-7r° background is about 20%. 

The contamination of the sample from sources of high energy single photons other than n° 
decays, such as prompt photon production, is expected to be negligible since the rate of such 
processes is low [44]. Using the Lepto and Ariadne models the total contamination due to 
the misidentification of electrons, and of ^-mesons decaying to two photons, is found to be less 
than 4%. The measurements are corrected for this using the QCD-based models following the 
procedure outlined in section 4.3. The contribution of the background from photoproduction 
processes was studied using the Phoiet [45] model and found to be negligible. 

4.3 Correction Procedure and Systematic Uncertainties 

The results shown in this paper consist of two sets of spectra. First, the dependence of the 
ep cross section for inclusive forward 7r°-meson production on Bjorken-rr, Q 2 , p* T lT and x n is 
studied. Measurements are then presented of the transverse energy flow in 7r°-tagged events. 
The transverse energy is evaluated from energy deposits measured in the LAr and SPACAL 
calorimeters, supplemented with tracking information from the CT, according to the prescrip- 
tion in [46]. 

The data are corrected using a bin-by-bin unfolding procedure using event samples gen- 
erated with the Lepto and Ariadne models. The correction factors are obtained by taking 
the average of the correction factors estimated by these two models. The typical values of the 
correction factors obtained are approximately equal to 1.5. 

The following sources of systematic uncertainty are considered. The uncertainties from 
each source are added quadratically to form the total point-to-point systematic uncertainties on 
each of the measured distributions presented here: 

• The model dependence of the bin-by-bin acceptance corrections leads to systematic un- 
certainties of between 4% and 11% on both the it cross section and transverse energy 
spectra. This source of uncertainty is calculated as half the difference of the correction 
factors derived from Lepto and Ariadne. The uncertainty resulting from the model 
dependence is largest at the lowest values of x and for the highest values of p* T lT and x n . 

• The systematic uncertainty due to the correction for QED radiative effects, calculated 
using Heracles, is typically 3% for both the it cross section and the transverse energy 
flow measurements. 

• Variation of the n° selection and identification cuts within the resolution of the recon- 
structed quantities gives rise to an uncertainty of 5% to 10% in the measurements of both 
the cross section and the transverse energy flow. 

• The uncertainty on the electromagnetic energy scale of the SPACAL (1%) affects the 
reconstructed kinematics of the scattered positron. This results in uncertainties on the 
7r° cross section and the transverse energy flow measurements of typically 9% and 2%, 
respectively. 
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• The uncertainty on the electromagnetic energy scale of the LAr (3%) leads to an uncer- 
tainty on the 7r° cross section of 5% to 10%, but has a negligible impact on the transverse 
energy spectra. 

• The hadronic energy scale uncertainty on the LAr (4%) gives rise to an uncertainty on the 
transverse energy flow measurements of 4%, but has negligible impact on the ir° cross 
sections. 

• The uncertainty on the hadronic energy scale of the SPACAL (7%) results in uncertainties 
on the 7r° cross section and the transverse energy flow measurements of typically 2% and 
5%, respectively. 

• The uncertainty on the polar angle measurement of the scattered electron (0.5 mrad) has a 
small influence (below 2%) on the 7r° cross section and the transverse energy flow spectra. 

• The uncertainty on the determination of the trigger efficiency leads to a 5% uncertainty 
on the 7T° cross section measurements, but has a negligible effect on the transverse energy 
flow spectra. 

• The uncertainty on the luminosity measurement leads to a 1.5% uncertainty on the 7r° 
cross section measurement, but has no effect on the transverse energy distributions. 



5 Results 

5.1 Inclusive Forward 7r° Cross Sections 

The inclusive it cross section is measured differentially as a function of Q 2 and x, and as a 
function of p? n and x n , for 7r°-mesons produced in the range p^ > 2.5 GeV, 5° < 9 n < 25° 
and Xj, > 0.01. The DIS phase space is restricted to the kinematic range 2 < Q 2 < 70 GeV 2 
and 0.1 < y < 0.6. The cross section data presented in this section are also given in Tables 1 
and 2. 

The inclusive cross section da w / dx for p* T w > 2.5 GeV as a function of Bjorken x is shown 
in Fig. 2 for three intervals of Q 2 : 2 < Q 2 < 4.5 GeV 2 ; 4.5 < Q 2 < 15 GeV 2 and 15 < Q 2 < 
70 GeV 2 . The distributions rise with falling x except at values of x of around 10~ 4 in the lowest 
Q 2 region. This turnover is due to the limitations in the phase space imposed by the ir° and DIS 
event selection cuts. 

The predictions of five QCD-based models are compared with the data. Calculations from 
RAPGAP which implement DGLAP evolution for proton structure only, labelled DIR, fall sub- 
stantially below the data. The disagreement becomes more pronounced at lower values of 
Bjorken-rr. Calculations which assume virtual photon structure, marked DIR+RES, describe 
the data well, although it is necessary to use rather large renormalisation and factorisation scales, 
/i 2 = Q 2 + 4p|,, in order to get a sufficiently large resolved photon component. Using the same 
scale in RAPGAP a reasonable description of the azimuthal jet separation in a measurement 
of inclusive dijet production at low x in DIS [47] is obtained. Forward jet data are also well 
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described by DGLAP-based calculations which include virtual photon structure although dif- 
ferent measurements tend to prefer different choices of scales [3,4, 26], fi 2 = Q 2 + \.p\ or 
fx 2 = Q 2 + p\. Using the latter value would result in a reduction of about 30% in the normal- 
isation of the predicted forward ir° distributions measured in this paper, with little change to 
the shape. Predictions based on the CCFM equation, labelled CCFM (CASCADE), agree at 
the highest values of Q 2 . However, they fall below the data at low values of x and Q 2 . The 
rate of forward jet events predicted by CASCADE agrees with the data [26]. This observation 
is not in contradiction with the forward n° data since the discrepancies observed here arise at 
the lowest values of x, which are not covered by the forward jet measurement. Differences 
between Rapgap and Cascade in the overall description of forward jet and particle produc- 
tion may be related to differences in the modelling of the partonic structure of forward jets. 
Cascade mostly produces gluon-induced jets, while Rapgap has a substantial contribution 
of quark-induced jets in the forward region. Fewer high momentum particles are produced in 
gluon-induced jets than in quark-induced jets [48]. Predictions of the CDM give a reasonable 
description of the data. Analytical calculations using a modified BFKL equation, labelled mod. 
LO BFKL, describe the data well in the lower Q 2 region although they have a tendency to 
exceed the data at the lowest values of Bjorken-rr in the highest Q 2 intervals. 

Fig. 3 shows the differential cross section da w /dx for transverse momenta p* T n > 3.5 GeV, 
in three intervals of Q 2 : 2 < Q 2 < 8 GeV 2 ; 8 < Q 2 < 20 GeV 2 and 20 < Q 2 < 70 GeV 2 . 
Compared with the spectra in Fig. 2, the measured differential cross sections are lower by factors 
of between two and four. The QCD-based models provide a similar quality of description to 
that given in Fig. 2. Calculations based on the modified BFKL equation clearly exceed the data 
in the two highest Q 2 intervals. The prediction of the NLO calculation [34] is also shown and 
describes the data well. Predictions of the CDM (not shown) overestimate the data at the lowest 
values of Q 2 and x. 

The cross section da n / dp^ w is shown as a function of p^ )7r in Fig. 4 in the same Q 2 intervals 
as in Fig. 2. The data fall steeply with increasing p* Tix and the shapes of the distributions 
vary only slightly with increasing Q 2 . The calculations implementing resolved virtual photons 
describe the data well. The predictions of the model including only direct processes fall below 
the data everywhere although they come nearer to the data as Q 2 increases. The CCFM-based 
calculations fail in the lowest Q 2 interval but give a reasonable description of the highest Q 2 
interval. The CDM predicts spectra which are somewhat harder than the data. 

In Fig. 5 the cross section da^/dx^ is shown as a function of x w in the same three intervals 
of Q 2 as used in Fig. 2. In Fig. 6 the differential cross section is presented in three intervals of 
Bjorken-rr: 4.2 • 10~ 5 < x < 2 • 10~ 4 ; 2 • 10" 4 < x < 10~ 3 and 10~ 3 < x < 6.3 • 10" 3 . The cross 
section falls as x w increases. There is no strong dependence of the shapes of the distributions on 
Bjorken-rr or Q 2 . The resolved photon approach, the CDM, and the BFKL calculations describe 
the spectra well. The CCFM implementation describes the data only in the highest intervals of 
Q 2 and x. The direct photon calculations fall below the data in all of the spectra but approach 
them as Q 2 increases. 
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5.2 Transverse Energy Flow 



The transverse energy flow, defined as the mean transverse energy per event per unit of pseu- 
dorapidity difference in the hadronic centre of mass system, -^dE^/d(rj* — rj*), in events con- 
taining at least one forward ir° is presented in Fig. 7. Here N is the total number of events and 
is the sum of the transverse energies of each particle i: E^ = H^E^. The transverse energy 
E^i of a particle i with energy E* and polar angle 9* is defined as E^ = E* sin 9*. The pseu- 
dorapidity r/* is defined as — lntan(#*/2). The energy flow, which includes the contribution 
from the forward ir°, is plotted as a function of the difference in pseudorapidity rj* — r]* from 
the selected forward 7r°. In events containing more than one forward ir°, the candidate with the 
largest transverse momentum is chosen. The spectra are presented in three intervals of the ir° 
pseudorapidity ranging from close (Fig. 7a) to far (Fig. 7c) from the proton direction. 

The spectra show a large increase of transverse energy production in the region associated 
with the 7r°. This can be understood as being due to the energy associated with the jet which 
contains the leading ir°. A broad distribution of lower transverse energy flow in the current 
region reveals the range over which the transverse momentum of the jet is compensated. 

The QCD-based models all describe the transverse energy flow in the vicinity of the ir° but 
give different predictions in the current region. Calculations with resolved photon processes 
tend to agree best with the data. The CCFM approach provides a reasonable description of the 
data. This model predicts a strong compensation of the ir° transverse momentum in the rapidity 
region between the forward particle and the proton remnant system which is, however, not 
covered in the measurements presented here. The direct photon model shows a peak at larger 
values of pseudorapidity difference rj* — r?* than is observed in the data. This effect becomes 
less pronounced with increasing ir° pseudorapidity as the forward it may enter the current jet 
region. The differences between the models can be qualitatively understood as a consequence 
of the ordering or otherwise of the kx in the parton cascades. The k T ordering in the direct 
photon model forces the compensation of the transverse energy in the current region while 
models without this requirement allow compensation close to the 7r°-meson. The CDM (not 
shown) predicts too much transverse energy in the vicinity of the ix° . This is a consequence of 
the overly hard ir° transverse momentum distribution predicted by the CDM and shown earlier 
in Fig. 4. 

Transverse energy flow in the region away from the forward n° and into the current region is 
further studied as shown in Fig. 8. The mean transverse energy over the region 1.0 < rj* — r]* < 
3.0 is plotted as a function of Bjorken-rr for the same three intervals of the ir° pseudorapidity as 
in the previous figure. The data show no significant dependence on Bjorken-rr although there is 
a tendency for the mean transverse energy to fall as pseudorapidity of n° approaches the current 
region. With the exception of the direct photon model, all of the QCD-based models give a 
reasonable description of the data. The direct photon prescription is only able to describe the 
data in the pseudorapidity region closest to the proton remnant. 

6 Summary 

Measurements are presented of ep interactions containing high transverse momentum forward 
going 7r°-mesons in the deep-inelastic scattering regime of 0.1 < y < 0.6, 2 < Q 2 < 70 GeV 2 
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and 4-10 5 < x < 6 ■ 10 3 . The inclusive ir° cross section is measured differentially as a 
function of Q 2 , Bjorken-rr and n° transverse momentum and energy for particles with pj, n > 2.5 
GeV, 5° < W < 25° and x n = > 0.01. The transverse energy flow relative to the direction 
of the forward ir° and the mean transverse energy in the vicinity of the forward ir° are also 
measured. 

The measurements presented are sensitive to the dynamics of parton evolution. Several 
different QCD-based approaches are confronted with the data. An approach implementing 
DGLAP evolution of proton structure underestimates the ir° cross section at low values of 
Bjorken-rr and Q 2 and overestimates the range at which the n° transverse momentum is compen- 
sated. Calculations implementing virtual photon structure provide the best description, albeit 
with a preferred choice of renormalisation and factorisation scale which is inconsistent with that 
required by other measurements of the hadronic final state in the forward region. Predictions 
based on CCFM evolution fail to describe the ir° cross section at the lowest values of Bjorken-x 
but give a fair description of the transverse energy distributions. The Colour Dipole Model gives 
a reasonable description of the n° cross section but predicts a transverse momentum distribution 
of 7r°-mesons which is significantly harder than is observed in the data. Calculations using next- 
to-leading order QCD matrix elements convoluted with ir° fragmentation functions describe the 
Bjorken-x dependence of the forward ir° cross section well. Predictions made using a modified 
BFKL equation describe the forward ir° cross section at the lowest Q 2 values but exceed the 
data at the highest Q 2 values. 
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Table 1: The inclusive 7r°-meson cross sections as shown in Figs. 2 and 3, presented together 
with the statistical and total uncertainties. 



P*T,n 


/ daTt \±stat / pb \ 




/ day \±stat / pb \ 


P*T,* 


/ da^ \±stat / pb \ 


2<Q 


2 < 4.5 GeV 2 


4.5 < Q 2 < 15 GeV 2 


15 < Q 2 < 70 GeV 2 


2.5—2.8 


21 ^ ±9 


2.5—2.8 


200^ 4 


2.5—2.9 




2.8—3.3 


146& 


2.8—3.4 


1 99 ±5 


2.9—3.5 


69.1||;| 


3.3—4.0 


60.0g 4 


3.4—4.1 


58^ 


3.5—4.7 


97 7±2.0 
z ' • ' ±4.7 


4.0—5.2 


21 7 ±L4 

zi - ' ±4.7 


4.1—5.2 


23.0±I;| 


4.7—8.0 


c 41 ±0.38 
°-^ i ±0.96 


5.2—8.0 


a o 7 ±0.38 
^•°'±0.83 


5.2—8.0 


r o±0.4 
°-°±1.2 


8.0—15.0 


n 9i ±0 - 04 

U - Zi ±0.06 


8.0—15.0 


n oq±0.06 
U - 2y ±0.26 


8.0—15.0 


n si ±0 - 06 

U - Oi ±0.12 







Table 2: The inclusive 7r°-meson cross sections as shown in Fig. 4, presented together with the 
statistical and total uncertainties. 
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Table 3: The inclusive 7r°-meson cross sections as shown in Figs. 5 and 6, presented together 
with the statistical and total uncertainties. 
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Figure 1: Generic diagrams for DIS processes at small x. The gluon longitudinal momentum 
fractions and transverse momenta are labelled Xi and kn, respectively, (a) A gluon ladder 
evolves between the quark box, attached to the virtual photon, and the proton, (b) The partonic 
structure of the photon is "resolved". 
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Figure 2: The inclusive ep cross section for forward it mesons produced in the range p* T-K > 2.5 
GeV, 5° < 9 W < 25° and x n = E n / E p >0.01 as a function of Bjorken-x in three intervals of 
Q 2 . The DIS kinematic region is further specified by 0.1 < y < 0.6. The inner error bars 
denote the statistical uncertainties and the outer error bars show the statistical and systematic 
uncertainties added quadratically. The predictions of five QCD-based models discussed in the 
text are shown. 
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Figure 3: The inclusive ep cross section for forward ir° mesons produced in the range p* T-K > 3.5 
GeV, 5° < 9 W < 25° and x n = E n / E p >0.01 as a function of Bjorken-x in three intervals of 
Q 2 . The DIS kinematic region is further specified by 0.1 < y < 0.6. The inner error bars 
denote the statistical uncertainties and the outer error bars show the statistical and systematic 
uncertainties added quadratically. The predictions of five QCD-based models discussed in the 
text are shown. 
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Figure 4: The inclusive ep cross section for forward n° mesons produced in the range p* T n > 2.5 
GeV, 5° < On < 25° and x n = E w / ' E p >0.01 as a function of pt,w in three intervals of Q 2 . The 
DIS kinematic region is further specified by 0.1 < y < 0.6. The inner error bars denote the 
statistical uncertainties and the outer error bars show the statistical and systematic uncertainties 
added quadratically. The predictions of four QCD-based models discussed in the text are shown. 
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Figure 5: The inclusive ep cross section for forward tt° mesons produced in the range -p* T n > 2.5 
GeV, 5° < 9 W < 25° as a function of x n in three intervals of Q 2 . The DIS kinematic region is 
further specified by 0.1 < y < 0.6. The inner error bars denote the statistical uncertainties and 
the outer error bars show the statistical and systematic uncertainties added quadratically. The 
predictions of five QCD-based models discussed in the text are shown. 
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Figure 6: The inclusive ep cross section for forward it mesons produced in the range p* T-K > 2.5 
GeV and 5° < 9 W < 25° as a function of x w in three intervals of Bjorken-a;. The DIS kinematic 
region is further specified by 2 < Q 2 < 70 GeV 2 and 0.1 < y < 0.6. The inner error bars 
denote the statistical uncertainties and the outer error bars show the statistical and systematic 
uncertainties added quadratically. The predictions of five QCD-based models discussed in the 
text are shown. 
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Figure 7: Distributions of transverse energy flow in events containing a forward ir° produced in 
the range p* T ^ > 2.5 GeV, 5° < 9 W < 25° and x n = E w /E p > 0.01. The transverse energy flow 
is presented as a function of the distance in pseudorapidity from the selected forward n° for 
various ranges in the ir° pseudorapidity. The DIS kinematic region is specified by 2 < Q 2 < 70 
GeV 2 and 0.1 < y < 0.6. The inner error bars denote the statistical uncertainties and the outer 
error bars show the statistical and systematic uncertainties added quadratically. The predictions 
of three QCD-based models discussed in the text are shown. 
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Figure 8: The mean transverse energy in events containing a forward ir° produced in the range 
p* T lT > 2.5 GeV, 5° < 9 n < 25° and x n = E n J ' E p > 0.01. The transverse energy is measured 
over the region 1.0 < rf — r/l < 3.0 as a function of Bjorken-a; for three intervals of tt° 
pseduorapidity. The DIS kinematic region is further specified by 2 < Q 2 < 70 GeV 2 and 
0.1 < y < 0.6. The inner error bars denote the statistical uncertainties and the outer error bars 
show the statistical and systematic uncertainties added quadratically. The predictions of four 
QCD-based models discussed in the text are shown. 
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